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Overview 

Nondestructive Evaluation Sciences Branch 

• Focus on composites 

• Validation related challenges 

• Experimental comparisons 

• Avoiding application of the tool to cases that don't fit prior validations 

• Tool development challenges 

• Memory efficiency 

• Computational efficiency (speed) 

• Hardware related challenges 

• Hardware selection and keeping up with continuous progress 



Composites for Space 
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Composites for Aeronautics 



• Advanced Composite Project (5 Year Project): 

• Reduce timeline for certification of composite structures 

• Partnership: NASA, FAA, DoD, Industry, University 

• Rapid Inspection Technical Challenge: 

• Focus areas: 
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• Inspection of complex geometry components 

• Rapid large area inspection 

• Damage/defect characterization 

• Validation of detectability 

• Damage types: 

• Microcracking, fiber waviness, delamination, porosity 

Simulation: 

• Enables model based inspection prediction/validation and 
cost effective method optimization 

• Custom code, 3D ultrasound simulation under 
development 


Composite Damage/Defect Types 



X-ray CT data of 
microcrack damage 



X-ray CT data of 
microcrack damage 



Micrograph showing resin rich 
regions and fiber misalignment 


Fiber waviness (in-plane). 

From Kugler and Moon 2002 
doi: 10.1177/0021998302036012575 


X-ray CT of PRSEUS Joint 
From NASATM-2013-217799 by Patrick Johnston 
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Micrograph showing porosity 


X-ray CT data of 
delamination damage 


UT data of 

delamination damage 
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Simulation 


Tool Validation Challenges 
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• Direct comparisons between simulation and experiment can be 
challenging 

• Requires specific experimental design 

• Experimental case will always have some differences from the simulated 
case 

• Getting representative samples for experiment can be a challenge 

• Creating representative defects/damage 

• Differences between 'idealized' material properties and as-manufactured 

• Must perform re-validation against appropriate cases when the simulation 
tool is used for a new purpose 

• Understanding of the physics is required to know when this is necessary 


Example: Ultrasound simulation 
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• Elastodynamic finite integration technique ultrasonic simulation code 
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• Custom C++ and MPI 

• Similar to finite difference 

• Adaptable, efficient, all details under our control 
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Simulation Validation Approach 1 
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• Laser Doppler vibrometry experiment comparisons 

• Group velocity comparisons unidirectional IM7/8552 8-ply sample: 

• Track envelope peak propagation (using Hilbert transform) 



a) LV wavefield data b) EFIT result 

Mode 1 group velocity comparisons: 


Method 

Mode 1 . 0° 
ty(m/s) 

Mode 1, 90° 
v g (m/s) 

% Difference 
from EFIT, 0° 

% Diff from 
EFIT, 90° 

EFIT 

1956 +/- 90 

1335 +/- 44 

- 

- 

Dispersion curves 

1911 

1254 

2.33 

6.26 

Experiment 

2254 +/- 84 

1464 +/- 69 

14.16 

9.22 


EFIT simulation result 
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Experimental LV wavefield data 



Note: time step of frames and colormapping in 
the two movies is not the same 






Simulation Validation Approach 2 

Nondestructive Evaluation Sciences Branch 

• Compare wave behavior in all directions 

• Wavenumber comparison technique: 

• Start with data for all grid points on surface of sample, amplitude at x-position vs. y-position 
vs. time 

• Take 3D FFT to yield x-wavenumber vs. y-wavenumber vs. frequency (where k=f/c phase ) 

• Select frequency slice that corresponds to the excitation frequency 
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Unidirectional (Ease 
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• AO mode results (SO mode is very low amplitude) 

• Amplitude variation around k x -k y oval is due to excitation source filtering (changes with couplant, transducer, etc and 
can be included in EFIT as well) 

• Interested in directional wave behavior observed via wavenumber values 



Wavenumber plots 



Simulation 



Unidirectional laminate: Mode 1 wavenumber comparisons 


Method 

0° jfe(l/m) 

90° k(l/m) 

% Difference 

% Diff from 




from EFIT, 0° 

EFIT. 90° 

EFIT 

143.3 

263.4 

- 

— 

Dispersion curves 

139.7 

258.1 


2.54 

2.03 

Experiment 

129.5 

229.1 


10.12 

13.93 


Cross-ply case 
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• Demonstrates ability to build up laminates ply-by-ply in EFIT 

• However - is orthotropic only! 
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Method 

0° k(l/m) 

90° k{l/m) 

% Difference 

% Diff from 




from EFIT, 0° 

EFIT, 90° 

EFIT 

159.8 

169.2 


- 

- 

Dispersion curves 

161.4 

180.7 


1.00 

6.57 


Experiment 

152.3 

170.5 


4.81 

0.77 




Orthotropic vs. Non-orthotropic 

Nondestructive Evaluation Sciences Branch 


rp(t+At/ 2) 

i ll 


rji(t — Z\t/2) 
= i 11 



|4?(4 n) - 4 n_il) ) + 4»(*4” > - 4“ _l2) ) + 4?(4 n) - 4"~***j] 



/ 2 ^ 1 2 ) 

-'ll —-'ll 


+ 


At 


,(w)( v (n) _ „,(n-£i) 


Hi 


-) + c<?( 4 ” ) 


Vi 


(n—x 2 )\ , (n) / ( n ) 




(n-xz)i 


+ 


4A.r 


,( n ) 

14 


(< 


„(«) _l_ „( n -*2) , ..(n+x 3) , (n-£ 2 +£3h 

v 2 ' v 2 ' v 2 "1 v 2 / 


(n—x 2 ) 
; 2 

(n-x 3 ) 


/ I7t— X2) , M \TI—X 2 —X 3 ) \n—x 3 J\ / 177 — X 31 mi . i 77 -t-X 2 — X 3 ; 1 n-\-X2)\ 

— (^2 + v 2 + v 2 + v 2 ) + \ V 3 + v 3 + v 3 + ^3 ) 


(71—X2—X3) 


(n— x 3 )> 


(n— £3) 


n 


( 77 +^ 2 )^ 


-(% 

+4? ((4 


. (n) ( 77 — X 2 ~ £ 3 ) , (n—x 2 ) 

+ V 3 + V 3 + ^3 


77 — Xl) 


+ Hi 71 '* + V 


(n— X1+X3) (n+x 3 ) 


1 


V 


1 


0 - K 


(n) . (n—x 1 — x 3 ) . (n— x 3 )> 

+ v\ ’ + Hi 31 + hJ 3; ^ 


+ K 


(n— £ 3 ) (n) (n+£i-£ 3 ) (n+£i) 

I ^3 I t /3 I t /3 


+>> 

' c 16 


((4 


+(4”~ i2> + 4”' + + 4"™') - M + <4 + v'2- + » 2 


(n— £ 1 ) (n— £ i+£ 2 ) , (n) (n+£ 2 ) 

v + hJ + hJ + H^ 


(n— £ 1 ) 

1 

(n— £ 1 — £ 3 ) (n-£i) 

I t /3 

(n-£i-£ 2 ) 


(n— £ 3 ) 


) - i v r ~ x + v f ~ i/ + v f + "3 


(n) 


,( n ) 


(n+x 1 —X 2 ) 


(n+x i)n 


) - K 

(n—x 1 — X 2 ) 


. (n— xi) (n) (n— X2) > 

+ Hi + Hj + Hi 


(n— £ 1 ) 
^2 


(n— £ 2 ) 

} 2 




Quasi-isotropic case 
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Differences from Experiment 







Overall simulation yields results close to experiment 
and predicted dispersion curves 

Expected some differences from experiment due to as- 
fabricated material properties of laminate versus "ideal" 
properties used in model 

• Thickness variation, fiber warping, variation in fiber density, 
slightly off-angle ply layers (laid up by hand) 


Unidirectional laminate: 


[0p/9(VL laminate: 



[0 2 /90 2 ] s laminate: 
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Once we have a working simulation tool - let's use it! 

Develop optimized and new damage quantification methods 
Predict inspectability 
Validate SHM 
Etc 

Still have validation challenges ahead 

You developed a promising new inspection methodology, now go 
back and validate against experiment 


Guided Wave Energy Trapping 



• Studied previously by several authors via LDV and simple simulations 
• These prior studies focused on single layer delamination 
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2 Sohn, H., Dutta, D., Yang, J. Y., Park, H. J., DeSimio, M., 
Olson, S., & Swenson, E. (2011). Composites science and 
technology , 71(9), 1250-1256. 



Glushkov, E, Glushkova, N, Golub, 
M, Moll, J, Fritzen, CP. Smart 
Materials and Structures 21.12 
(2012): 125001. 
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Michaels, Composite, simulated single delam 
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Michaels , J; Dawson, A ; 
Michaels, T ; Ruzzene, M. Proc. 
SPIE 9064, (2014); 
doi:10. 1117/12. 2045172. 
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Tian, Composite, single delam 



3 Zhenhua Tian ; Lingyu Yu ; Cara A. C. Leckey; 

Proc. SPIE 9063, (2014), doi:10.1117/12.2044927. 


Energy Trapping 



• Potential for rough sizing of damage via 


LDV data: 
500 kHz 


LDV data: 
200 kHz 



Mass normalized 
cumulative energy 

rh 1 

Ei(x,y,z, t) = I -vfdt 
h 1 A 



rapid data processing 
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• Can energy trapping be leveraged for multi-ply delam characterization 
with only single sided access? 





Hidden Delamination Study 



• Specifically: can energy reduction be used to 
provide information about the presence of 
hidden delamination damage? 

• Simulation based study: 

• 8 ply, IM7/8552 CFRP sample [(0/90) 2 ] s , 0.92 mm thick 

• 3 simple delamination cases: 1, 2, and 3 delaminations (+ pristine 
case) 

• 300 kHz, 3 cycle Hann windowed sine wave 

• dx=19 pm, dt analysis = 0.29 ps (dt/200) 
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Simulation Study: Cumulative Energy 
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Energy trends 



• Quantitative comparisons, energy above top delamination 



Note: nonzero for pristine 
because energy still passes 
into that region (especially 
with edge scattering, etc) 


v Sum of Normalized Cumulative Energy Above Uppermost Delamination 
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Validation 

• LDV scans of composites with Teflon inserts 

• Matching simulated cases (ranging from 
single delam to 3) 

• Contact transducer for excitation 

• Plan is to measure cumulative energy 
trapping 
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Inclusion with Teflon Case 
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• Unexpected defects/variability due to manufacturing 

• Immersion UT 10 MHz, 0.01" resolution 
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Repeatability 



• For amplitude based studies, repeatability can be an issue 

• Differences between repeated data scan sets 

• Couplant changes! 

• Also, planning to use laser excitation 


Transducer face: Scan 1 
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Transducer face: Scan 2 
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Code optimization 


• Single vs double precision? 

• Porting code often requires some 
amount of re-writes 

• Must check that ported code 
yields same results as validated 
code! 

• Example shown here: 

• ~10x improvement in 
efficiency 

• 0.6% cost 

• But - scalability also matters 
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Comparison of Platforms 


Optimized Code 



Original Code 


HPCl, 70 

~$150K 

8 nodes 
70 CPUs 


K2, 160 

~$800K 

32 nodes 
160 CPUs 


i7-IPP, 10 i7-lntrinsic, 10 3120A, 166 7120P, 166 7120P, 243 


~$2K 
1 node 
10 CPUs 


~$2K 
1 node 
10 CPUs 


~$3K 
1 node 
166 CPUs 


~$5K 

1 node 
166 CPUs 


~$5K 

1 node 
243 CPUs 


Conclusions 



• Inevitable shift towards more use of modeling and simulation to test 
hypotheses, optimize methods, predict inspectability, etc 

• This is a good thing, as it may enable new approaches and cost- 
effective investigation of inspection methods 

• However, validation is a key step with many challenges 

• Whether custom or commercial simulation software is used, it is of 
key importance to know that validation has been performed and 
successful 


Thanks to: Eric Burke, NASA LaRC 
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